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ABSTRACT 

We study the demographics of supermassive black holes (SMBHs) in the lo- 
cal and high-z universe. We use the continuity equation of the population of 
SMBHs as the leading principle. We consider three physical processes for the 
growth of SMBHs: mass accretion, mergers, and direct formation of SMBHs. 
The mass accretion history of SMBHs is estimated by the hard X-ray luminos- 
ity functions (HXLFs) of AGNs. First, we compare the mass accretion history 
at z > with optical luminosity functions (OLFs) of QSOs previously studied 
and that with HXLFs. We conclude that the constraints on parameters of mass 
accretion (energy conversion efficiency, e and Eddington ratio, /Edd) based on 
the continuity equation appear to be adequate using HXLFs rather than OLFs. 
The sub-Eddington case (/Edd < 1) is allowed only when we use HXLFs. Next, 
we extend the formulation and we can obtain the upper limit of the cumula- 
tive mass density of SMBHs at any redshifts. For an application, we examine 
if Lyman-Break galaxies (LBGs) at z ~ 3 already have SMBHs in their centers 
which is suggested by recent observations. We tentatively assume the presence 
of SMBHs in LBGs and that their mass, Mbh is proportional to the stellar mass 
of LBGs, M* with the mass ratio £ = Mbh/M*- If most of LBGs already has 
massive SMBHs at z ~ 3, the resultant mass density of SMBHs at z ~ should 
exceed the observational estimate because such SMBHs should further grow by 
accretion. Therefore, we can set the upper limit of the probability that one LBG 
has a SMBHs. Since the merger rates and direct formation rates of SMBHs are 
uncertain, we consider two limiting cases : (i) mergers and/or direct formations 
are not negligible compared with mass accretion and (ii) mass accretion is the 
dominant process to grow the SMBHs. The special conditions should be met in 
order that a large part of LBGs have SMBHs in both cases. In case (i), we may 
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assume the constant parameters of mass accretion of AGNs for simplicity. Then, 
large energy conversion efficiency and frequent mergers and/or direct formations 
at z > 3 are needed so that a large part of LBGs have SMBHs with £ = 0.002 to 
0.005. Whereas, in case (ii), energy conversion efficiency should be mass depen- 
dent and the constraint is strict; the fraction of LBGs which have SMBHs must 
be less than 10-40%. In both cases, the possibility that nearly all LBGs have 
SMBHs with large mass ratio, such as £ > 0.005, is reliably ruled out. 

Subject headings: black hole physics — galaxies: active — galaxies : evolution 
- galaxies: luminosity function, mass function — galaxies : nuclei — quasars: 
general 

1. Introduction 

Today, there is a wide consensus that almost all local galaxies harbor supermassive 
black holes (SMBHs) in their centers (e.g. Begelman 2003, Kormendy & Richstone 1995). 
We can estimate the mass of these SMBHs with dynamical motion of gases or stars around 
the center of galaxies (e.g. Ghez et al. 1998, Richsotne et al. 1998, Gebhardt et al.2003). 
It is well-known that there are some correlations between central SMBHs and its galactic 
bulges. For example, the mass of SMBHs is proportional to the mass of their galactic bulges; 
M B H/-^buige ~ 0.001 — 0.006, though its tightness and linearity is still under discussion 
(Kormendy & Richstone 1995, Maggorian et al. 1998, Laor 2001, McLure & Dunlop 2002, 
Marconi & Hunt 2003). Another tight correlation is so-called M B h — o relation, where a is 
the stellar velocity dispersion in the bulges. This relation is expressed as the power-law form, 
M BH oc a 3 - 8 - 4 - 8 (e.g. Gebhardt et al. 2000, Ferrarese & Merritt 2000, Tremaine et al. 2002). 
Recently, Ferrarese (2002) and Baes et al. (2003) argue that SMBH mass is related to the 
dark halo mass beyond the bulge, Mbh oc M^ 7 " 1 ' 82 . Graham et al. (2001) point that there 
is another correlation that bulges with more massive SMBHs have steeper central cusps. In 
spite of discoveries of many correlations and great efforts which have been devoted to clarify 
the physical meanings of these correlations (e.g. Silk & Rees 1998, Ostriker 2000, Umemura 
2001, Adams, Graff & Richstone 2001), we do not have a consensus yet regarding physical 
link between the formation process of SMBHs and that of their host galaxies (bulges). 

To understand how the central SMBHs formed and how their formation history is related 
to that of galaxies, we must at least know when SMBHs were formed in galaxies. Therefore, 
it is very important to know whether high- z (not active) galaxies already had SMBHs or 
not. In this respect, it is interesting to note that the so-called Lyman-Break technique, for 
example, has made it possible to discover a significant number of high- z (> 3) galaxies, i.e., 
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Lyman-Break Galaxies (LBGs, Steidel et al. 1996). These LBGs are widely known as high-z 
starburst galaxies and their typical star for mat ion rate is 10 2 — lO 3 M yr~ 1 . Papovich et al. 
(2001) and Shapley et al. (2001) applied the population synthesis models to infer UV-to- 
optical spectrum of LBGs at their rest frame and examined various properties of LBGs. 
Despite uncertainties their results are consistent at some points; LBGs are typically young 
(their ages are several hundred Myr since last star-formation events) and small (their typical 
stellar mass is ~ 10 10 M), though there were some old LBGs. 

Ridgway et al. (2001) point that the sizes and magnitudes of the host galaxies of the 
radio-quiet QSOs at z ~ 2 — 3 are similar to those of LBGs in HST imaging study. Steidel 
et al.(2002) observed about 1000 LBGs and argue that 3% of these LBGs are optically 
faint AGNs based on their UV spectrum. Nandra et al.(2002) confirmed by using Chandra 
observatory that the same number fraction, 3%, of LBGs are extremely bright in the rest- 
frame hard X-ray band, and asserted that these LBGs show AGN activity. They discussed 
that this number fraction, 3% may reflect the duty cycle of mass accretion to SMBHs in 
LBGs. That is, if the typical age of LBGs is 300 Myr (e.g. Shapley et al. 2001) and 
the accretion timescale is 10 7 yr, remaining 97% of LBGs may have inactive SMBHs. This 
situation may be compared with the local universe. In the local universe, almost all galaxies 
harbor SMBHs in their centers, but there are much inactive galaxies than AGNs. The 
pioneering work (Dressier, Thompson & Shectman 1985) found that the number franction 
of AGNs to all galaxies is only several %, though this franction recently increases to 20-40 % 
due to the accurate observations (Kauffmann et al. 2003, Miller et al. 2003). Furthermore, 
Vestergaard (2003) and Vestergaard (2004) discusses that the LBGs with AGN activity may 
have SMBHs of order 1O 8 M from the virial mass estimate using C IV line width. If these 
LBGs have typical stellar mass of M* ~ 1O 1O M , the mass ratio becomes Mbr/M* ~ 0.01, 
which is larger than MBH/^buige ~ 0.001 — 0.006 of local galaxies. Incidentally, normal 
galaxies and AGNs have the same mass ratio, MeH/^buige, i n the local universe (Wandel 
2002). Whereas, Granato et al. (2001) found that LBGs may be equivalent to the galaxies 
which are in the pre-AGN phase based on their semi-analytic model of the joint formation 
of QSOs and galaxies. Kawakatu et al.(2003) argue that it may be physically difficult to 
feed a black hole in optically thin galaxies, such as LBGs, based on the idea by Umemura 
(2001). In the local universe, almost all the SMBHs observationally accompany the stellar 
spheroidal (bulge) component (M33, which is the bulgeless galaxy, appears to have no SMBHs 
(Gebhardt et al. 2001)), therefore, the bulge component may be indispensable to form a 
SMBH. If LBGs are still in the initial phase of forming the bulge component (Matteucci & 
Pipino 2002), LBGs may not have SMBHs. 

Of course, it is difficult to resolve the central parsec scale region of LBGs and search 
the (inactive) SMBHs observationally. In this paper, we theoretically examine the possi- 
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bility that LBGs have SMBHs from the viewpoints of the demographics. We consider the 
cosmo logical accretion history derived from the observed hard X-ray (2-10keV) luminosity 
functions (HXLFs) of AGNs at z > (Ueda et al. 2003, hereafter UAOM03). Yu & Tremaine 
(2002) (hereafter YT02) use the continuity equation for the number density of SMBHs (see 
also Soltan 1982, Small & Blandford 1992) and show that the accreted mass estimated with 
optical LFs (OLFs) of QSOs at z > is comparable with the SMBH mass density of the 
local universe. However, we can see only a part of various populations of AGNs in the optical 
band, thus the analysis with hard X-ray LFs, which include absorbed AGNs, is important. 
We consider the cosmological accretion history estimated with HXLFs and compare it with 
the accretion history estimated with OLFs (YT02). Furthermore, we basically extend the 
formulation of YT02 and crudely estimate how many and how massive SMBHs are allowed 
to exist in the high- z universe. If most LBGs already had very massive SMBHs at z ~ 3, 
the accreted mass into SMBHs until z ~ should exceed the BH mass density of the local 
universe. In this way, we will be able to set constraints on the mass density of SMBHs in 
LBGs. 

Below, we use the continuity equation of SMBH population as the leading principle. In 
§2, we first transform the continuity equation to the useful integrated form (§2.1). We then 
compare the cosmological accretion history at z > with the local population of SMBHs 
and constrain the physical parameters; the energy conversion efficiency, e and the Eddington 
ratio, /lbg(= -^boi/^Edd) (§2.2). In this time, we compare the accretion history with HXLFs 
and that with OLFs. In §3, we extend YT02's formulation and derive the upper limit of the 
population of SMBHs at z ~ 3 (§3.1). Furthermore, we consider the tentative SMBHs in 
LBGs with mass ratio, £ = Mbh/M* (M» is the stellar mass of LBGs), and crudely estimate 
the upper limit of the probability that one LBG has a SMBH (§3.2). In §4, we summarize 
and discuss the results. Throughout this paper, we adopt the A-dominant universe and the 
cosmological parameters are (Q M ,Q\,h) = (0.3,0.7,0.7). 

2. Cosmological Accretion History v.s. Local Population of Supermassive 

Black Holes 

2.1. Continuity Equation of the Population of SMBHs 

We begin with the continuity equation for the number density of SMBHs: n(M,t) = 
dN/dM(M,t); that is 

dn d[n- (M)] , , 

7i7 + 7th = 7me rg c(M, f) + 7form(M, t) (1) 
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(Soltan 1982, Small & Blandford 1992, YT02). Hereafter, we simply express the black hole 
mass as M. The source terms on the right-hand side, 7mer g e(-^M) an d 7form(^, t), represent 
the SMBH-SMBH merger rates (e.g. Haehnelt 1994) and the rates of the direct formation 
of SMBHs by, say, collapse of the supermassive stars (SMS) (~ 1O 5 M , or may be more 
massive). We can express these terms explicitly. The merger term is, 

Tmerge (-^"j ^) 

-ldM l ldM 2 n(M 1 ,M 2 ) [5(M - Mi - M 2 ) - 5(M - Mi) - 5(M - M 2 )] , (2) 
2 Jo Jo 

where tt(M 1 , M 2 ) is the distribution function which represents the merger rates among the 
SMBHs of mass Mi and M 2 . The factor 1/2 is necessary to correct the double counts of 
the merger rates in the integral in eq.(2). Several recent studies investigate the distribution 
function of mergers, f2(Mi,M 2 ) with the semi-analytic model of galaxies and QSOs (e.g. 
Kauffmann & Haehnelt 2000, Volonteri, Haardt & Madau 2003). However, the key problem 
of the SMBH merger is the timescale. That is, the timescale to refill the "loss cone" of 
the SMBH binary (see Begelman et al. 1980 or Milosavijevic & Merritt 2003). Although 
the hardening timescale of the binary have been actively studied, this problem is still under 
discussion (Makino 1997, Milosavijevic & Merritt 2001, Milosavijevic & Merritt 2003, Makino 
& Funato 2004). Therefore, we treat Q(M±, M 2 ) as the unknown function. Here, we assume 
that the total mass of SMBHs is conserved through merger events. As for 7f orm , we only 
note that this numerical value is not negative, 

7form (M,t) > 0. (3) 

Recently, the gravitational collapse of SMS is being studied with the general relativistic 
numerical simulations (e.g. Shapiro & Shibata 2002, Shibata & Shapiro 2002). However, it 
is still uncertain how and when the initial condition of the simulations is realized. Eventually, 
these source terms are highly uncertain and we do not know whether these terms exist or 
not (the merger term will be at least non- negligible (e.g. Merritt & Ekers 2003)). YT02 
successfully avoid this difficulty by multiplying eq.(l) by the function, 

f(MM') = { °' forM '< M (4 ) 

J {M,™ ) | M'-M, for M' > M { ] 

and integrate by M' from to oo and by t from to t (cosmological time at z — 0), assuming 
that n(M, t) — > in the limit of t — > and M' — > oo. The result is 

Glocal(^f, to) = G acc (M, t ) + Gmcrge+form^, ^o)- (5) 

Here, G ioca i(M, t ), G acc (M,t ) and G mcTgc+iorm (M,to) are defined by 

POO 

Giocai(M,t )= / (M'-M)n(M',to)dM', (6) 

J M 
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G acc {M, t ) = / dM' I ° n{M', t)(M') dt, (7) 



M JO 



roo rt 

^mergc+form 

(M, t ) = / / (M' - M) (7mcrgc(M', t) + 7form(M', t)) dM'dt, (8) 

Jm Jo 

where (M') is the average mass accretion rate into the SMBHs of mass M'. Since G , mcrgc+ f or m(^, to) 
is or positive, we obtain the inequality, 

Gi ocal (M,t ) > G acc {M,t ). (9) 

These mathematical derivation of eq.(9) is the same as that in YT02, whereas we explain 
the physical meaning of eq.(9) here. To do so, we transform inequality (9) as 

poo poo rto poo 

/ M'n(M',t ) dM' > dM' n(M' ,t){M')dt + M n(M',t )dM'. (10) 
Jm Jm Jo Jm 

In this form, we can understand the meanings of each term. Clearly, the left-hand side of 
eq.(10) represents the mass density contained by SMBHs more massive than M at z = 0. 
The first term of the right-hand side represents the mass accreted into the SMBHs more 
massive than M per unit volume during < t < to (i-e. z > 0). The second term on the 
right-hand side represents another contribution; that is, the mass density of SMBHs at the 
moment that SMBHs become more massive than M. In our consideration, the SMBHs more 
massive than M can form through the mass accretion, mergers, or gravitational collapse 
of the SMS. The mass of a SMBH continuously grows by the mass accretion, whereas the 
mass of a SMBH grows discontinuously by the mergers and/or direct formation. Therefore, 
the mass of SMBHs at the moment that SMBHs become more massive than M is equal to 
M for the case with the mass accretion, whereas it is more than M otherwise. The total 
number of SMBHs which became more massive than M until z = is necessarily larger than 
n(M', to)dM' since the mergers reduce the number of SMBHs. Totally, the lower limit 
of the contribution of the mass density of SMBHs at the moment that SMBHs become more 
massive than M is M n(M' ,t )dM' ', which is the second term on the right-hand side of 
eq.(10). The physical meaning of this term may be considered as the "advection" of SMBH 
mass across the mass M in terms of the continuity equation of the population of SMBHs. 
This "advection" term is the heart of this inequality. 

We can calculate Giocai(^, ^o) with the local mass function of SMBHs. We adopt the 
mass function of SMBHs given by Aller & Richstone (2002) (hereafter AR02) based on the 
luminosity functions of galaxies at z = for each morphology. They convert the local LFs of 
normal galaxies to the mass function of SMBHs with some empirical relations among total 
luminosity, bulge luminosity, velocity dispersion and SMBH mass. One of these relations 
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is M — a relation mentioned above. AR02 derive the two local mass functions of SMBHs 
for different galactic LFs ; mass function (a) is derived from LFs by Marzke et al. (1994), 
whereas mass function (b) is derived from LFs by Madgwick et al. (2002). Below, we show 
the results of the calculation for the cases with mass function (a) and refer the results of the 
other cases in text. We use the fitting formulae for the mass functions of SMBHs given by 
AR02. With the mass functions (a) and (b), the local mass density of SMBHs more massive 
than 1O 6 M is (4.8 ± 1.6)h 2 x 10 5 M o Mpc~ 3 and (6.9 ± lA)h 2 x lO 5 M Mpc" 3 respectively. 
YT02 adopt the LFs based on the SDSS (Sloan Digital Sky Survey) and their calculated 
value is (2.5 ± 0.4) x 10 5 M o Mpc~ 3 with h = 0.65. Salucci et al. (1999) and Merritt & 
Ferrarese (2001) also obtain similar results. 

To calculate G acc (M, t Q ), there is a difficulty in that we do not know the number density 
of SMBHs at cosmological time t, n(M, t), and the average accretion rate of SMBHs, (M). 
Here, we assume that the mass accretion into SMBHs occurs mainly at bright AGN phase. 
The average accretion rate of AGNs is assumed to be 

(M) = ii^il^, (11) 

where L^ \ is the bolometric luminosity of AGNs and e is the energy conversion efficiency 
With eq.(ll), we can express G acc (M,t ) as 



cwM.^r ^r ^(f;"(i-f M ,, (12) 

Jl(m) Jo dU ec 2 



Here, d®(L,t)/dL is the LFs of AGNs and we use observed HXLFs (UAOM03). Very 
recently, UAOM03 build the HXLFs of AGNs over the wide luminosity range, \og(L x (2 — 
lOkeV) erg/s) = 41.5 — 46.5, and wide redshift range, z = — 3, from the combination of 
surveys with HEAOl, ASCA and Chandra satellites. The shape of HXLFs of AGNs at z = 
are fitted by so-called double power-law function well. That is, 

dH , L r\ = Q) = A l^/^r' 1 + ( w^*H - 1 • (is) 

d log L x 

In the above expression, 7I and 72 represent the inclinations of LFs at the lower and higher 
luminosity ends respectively, and represent the luminosity of the "knee" of LFs. The 
redshift evolution of HXLFs are described as the luminosity-dependent density evolution 
(LDDE), which is suggested by Miyaji et al. (2001) for the evolution of the soft X-ray 
luminosity functions of AGNs. The LDDE model is expressed as 



(14) 
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where 

, T v J (l + ^) pl , for * < ^ C (L X ) n , 

elZ '^ J 1 e(z c )[(l + z)/(l + z c (L x ))p> 2 , forz>z c (L x ), UDj 

and 

/ r n [ z* c , for L x > na\ 

^ = { K (L X /L a r, forix<L„. <16) 

The parameters which represent redshift dependence of e(z, L x ) are pi > and p2 < 0. 
Then, z c (Lx) corresponds to the peak redshift at which the number density of AGNs with 
luminosity of L x is the maximum. The cut-off redshift, z c (Lx), increases as the luminosity 
increases (a > 0) (see Fig. 10 of UAOM03). That is, the number density of the more luminous 
AGNs reaches the maximum at the higher redshift and decrease as z — > (see Fig. 12 of 
UAOM03). We simply extrapolate the fitting formulae given by above LDDE model at 
z > 3.0. 

To calculate G acc (M,t ), we need the relation between the hard X-ray luminosity, L x 
and the bolometric luminosity, L ho \. We assume the relation as L ho \ = CxLx, where C x 
is constant value C x = 25.0 (Alonso-Herrero et al. 2002). Generally, the ratio of B-band 
(4400 A), hard X-ray (2-10 keV) luminosity to the bolometric luminosity, C-q = L-q/L\ m \ 
and Cx should depend on the mass accretion rate and the black hole mass. This may be 
important to compare the local mass function of SMBHs with the mass accretion history 
estimated with OLFs and that estimated with HXLFs. We will discuss this issue in the 
next subsection in some details. We suppose that the QSOs radiate as L bol = /Edd-^Edd = 
1.5 x 10 38 /Edd(^A^©)(erg/s). The parameters here are the energy conversion efficiency, e, 
and the Eddington ratio, /Edd- 



2.2. Comparison between the Mass Accretion History with HXLFs and that 

with OLFs 

Here, we compare the mass accretion history estimated with HXLFs and that with OLFs 
based on the continuity equation. For both cosmological accretion history with HXLFs and 
that with OLFs, the relation to the local mass function of SMBHs can be investigated with 
ineq.(9) derived from the continuity equation. 

YT02 noted that their calculated G , acc (M,t ) with OLFs of QSOs with the typical 
parameter, e = 0.1 and /sdd = 1-0 is larger than G\ oca \{M, to), especially in the large mass 
range, M > 1O 7 M (inequality (9) is not satisfied). That is, there is the problem of the 
overaccretion compared to the local mass function of SMBHs. YT02 discussed the possible 
origin of this discrepancy. For example, they argue that the optically bright QSOs may 
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radiate at large mass-to-energy conversion efficiency, e (e.g. e > 0.2). However, we can 
see only a part of the light emitted in the mass accretion process to SMBHs in the optical 
band. This is partly because the optical luminosity emitted through the mass accretion 
can be absorbed by obscuring torus (Type II AGNs) (e.g. Antonucci 1993). We can miss 
the reddened QSOs due to the dust extinction in the host galaxies when the typical color 
selection criteria (e.g. UV excess) is adopted to build the QSO LFs (e.g. Webster et al. 
1995). Thus, we calculate G acc (M, t ) with HXLFs and compare it to G acc (M, t ) with OLFs 
and to G\ oca ,\(M,t ). Of course, as YT02 noted, it seems to make the problem worse to use 
the HXLFs. This is because HXLFs include extra populations of AGNs besides the optically 
bright AGNs. However, this is not the case. The redshift evolution of HXLFs and that of 
OLFs are different, then we can make the situation better (see below). 

The upper panel of Fig.l represents the results with the local mass function of SMBHs 
(a). For OLFs of QSOs, we adopt ones given by Boyle et al.(2000) at 0.35 < z < 2.3 derived 
from about 6,000 QSOs detected by 2QZ survey and simply extrapolate it at z > 2.3. We 
assume the constant bolometric correction, C B = 11.8 (Elvis et al. 1994). YT02 use the 
same OLFs and adopt the same assumptions. Like YT02's result, our results show that the 
inequality (9) is not satisfied in all mass range (especially at high mass end) with OLFs 
and typical parameters, (e, /Edd) — (0.1,1.0). If we use mass function (b), the calculated 
G\ oca ,\(M, to) is larger than (a) by a factor 1.4, but equation (9) is not still satisfied at 
M BH > 1O 7 - 5 M . 

However, Fig.l also shows that calculated G acc (M, t ) with HXLFs with the same pa- 
rameter set do not exceed G\ oca \(M,to) in M < 1O 8,5 M . At first sight, this looks strange. 
Since in HXLFs optically absorbed AGNs and reddened AGNs are considered besides the 
optically bright QSOs, G acc (M, t ) calculated based on HXLFs should become larger than 
G acc (M, t ) based on OLFs. The origin of this puzzle lies in our assumption on C B and 
C x - We assume that both of C B and C x are constant; C B = 11.8 and C x = 25.0. If this 
were actually the case, the ratio between B-band luminosity and hard X-ray luminosity, 
L-q/Lx, must be constant. However, UAOM03 compare their HXLFs of only optical type-I 
AGNs to OLFs by Boyle et al. (2000) and find that Lx is not proportional to L B . They 
obtained a power-law relation between the 2keV and 2500A luminosities, lx oc / Q 7 in the 
(£l M ,£l A ,h) = (1.0,0.0,0.5) universe. This relation is approximately transformed to the 
relation between L x (2 — lOkeV) and M B as 



That is, L B oc L^ 4 . Therefore, the assumption that both of C B and Cx are constant is not 
justified. As noted by UAOM03, if L B /Lx were constant, the redshift evolution of OLFs 
would be more quick than that of HXLFs. The bright end of OLFs grows more quickly than 




(17) 
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HXLFs as z increases. Then, calculated G acc (M, to) based on OLFs is larger than that based 
on HXLFs. Kawaguchi, Shimura & Mineghige (2001) build a disk-corona model for the 
broad-band spectral energy distributions (SEDs) of AGNs. Based on their model, Hosokawa 
et al. (2001) and Alonso-Herrero et al. (2002) calculated how Cb and Cx depend on the 
black hole mass, M and on the mass accretion rate, M. Though their dependence on M 
and M is not simple, Cx comparatively remains constant for varying M and M. Therefore, 
we adopt constant C x = 25.0 as noted above and use this value for the calculations below. 
Incidentally, if L x is proportional to the bolometric luminosity, L ho \, eq.(17) may be expressed 
as Cb oc L h ^ 4 oc M~ 0A with constant Eddington ratio. Alonso-Herrero et al. (2002) show 
the similar dependence, C B oc M" a25 based on Kawaguchi, Shimura & Mineshige (2001), 
though the dependence is somewhat weaker. 

Even if we use HXLFs, the calculated G acc (M, t ) with (e, / E dd) = (0.1,1.0) becomes 
larger than G\ oca \(M, t ) in the high mass end (M > 1O 8 ' 5 M ). However, this may be an 
inevitable result. We calculate G\ oca _i(M,t ) with the local mass function of SMBHs derived 
with the Schechter-type LFs of galaxies, which declines exponentially at the bright end, 
whereas we calculate G acc (M, to) with the HXLFs which is power-law function at the bright 
end. YT02 point that the QSO LFs at the bright end or M — a relation at the high mass 
end may be uncertain. Then, we use the energy conversion efficiency, e, and Eddington 
ratio, /lbg as constant free parameters to calculate G acc (M, to), and search the region of the 
parameter space of (e, /lbg) where inequality eq.(9) is satisfied only at M < M for given 
M . We choose M = 1O 8 -°M and M = 1O 8 - 5 M . Of course, the allowed parameter space 
becomes narrower as M increases, since the discrepancy we show in Fig.l is more prominent 
at the high mass end. 

Fig. 2 represents the results for the cases with the local mass function of SMBHs (a). 
The allowed region for each M is above the labeled lines respectively. We set the upper 
limit of e as e = 0.3, which is the maximum value for the thin-disk accretion models (Thorne 
1974). If we use mass function (b), the allowed region for each Mq becomes a little wider 
than (a). This is consistent with that the calculated G\ oca i(M,to) with (b) is larger and it is 
easier to satisfy the inequality (9) than the cases with mass function (a). As shown in Fig.2, 
to satisfy ineq. (9) at M < 1O 8 ' 5 M , the allowed region is not so large if we demand f^dd < 1- 
If /Edd ~ 1, e can be 0.1-0.3. However, only the large energy conversion efficiency (e ~ 0.3) 
is possible as /Edd decreases (/Edd ~ 0.3). This conclusion will not change, even if we adopt 
mass function (b). In Fig.2, we show the same parameter regions when we use OLFs of QSOs 
to calculate G acc (M, t ). To satisfy inequality (9) to high mass end, high Eddington ratio 
(/Edd > 1) and high energy conversion efficiency (e > 1) are needed. As shown above, the 
cosmological accretion history estimated with HXLFs naturally satisfy ineq. (9) to high mass 
end with typical parameters (especially /sdd < 1) and constant luminosity ratio, Cx = 25.0. 
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Therefore, we conclude that the accretion history with HXLFs is more plausible than that 
with OLFs. Below, we continue to use HXLFs to estimate the cosmological accretion history 



In the lower panel of Fig.l, we plot the calculated G acc (M, t ) at points 2 and 3 in Fig.2. 
As shown in Fig.l, if we adopt constant parameters, the quantity, G acc (M, to) — G\ oca \(M, to), 
increases as Mbh decreases to satisfy ineq. (9) upto high- mass end (large M ). From eq.(5), 
this means that the contribution of mergers and/or direct formation of SMBHs becomes 
important at the low mass end. If more massive halos possess more massive SMBHs, this 
may be the case since mergers between small halos occur more frequently compared with 
mergers between large halos under the CDM cosmology. 

In the next section, we consider the possibility that LBGs have inactive SMBHs as in the 
case of the normal galaxies at z — 0. To do so, we calculate the accreted mass into SMBHs 
during < z < 3. In this case, we should choose parameter set of (e, /lbg) so that ineq. (9) 
is satisfied. Therefore, we use (e, /lbg) indicated by points in Fig.2. As discussed above, 
with constant parameters, contribution of the mergers and/or direct formation of SMBHs 
to grow SMBHs becomes important in the low mass end. Below, we adopt the constant 
parameters for the case of non-negligible contribution of mergers and/or direct formation of 
SMBHs compared with mass accretion. 



3.1. Upper Limit of the Cumulative Mass Density of SMBHs at High-;? 



In this section, we investigate the possibility that LBGs have SMBHs. Our approach is 
the statistical one. That is, if there were already too many or too massive SMBHs in the 
universe of z ~ 3, the resultant mass density of SMBHs at z ~ should exceed the observed 
local mass density of SMBHs, that is a contradiction, because SMBH mass density can never 
decrease. The basic formulation is similar to that presented in the previous section. Here, 
we multiply the continuity equation, eq.(l), by f(M,M') defined by eq.(4) and integrate it 
by M' from M to oo and by t from t z=3 (cosmological time at z — 3) to to. We obtain 



of SMBHs. 



3. Demographics of SMBHs in the High-^ Universe 



Universe 



/ 

J a 



M 



oo 



M' (n(M', t ) - n(M', t z=3 )) dM' 




M— dM'dt 
dt 
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poo pto 

+ / / (M'-M) ( lmerge (M\ t) + 7 form(M', t)) dM'dt, (18) 
Jm Jt z=3 

where n(M, t ) is the local mass function of SMBHs, dN/dM(M,to), as we defined in the 
previous section. We substitute eq.(l) for the second term on the left-hand side of eq.(18), 
we obtain the equation, 

F local (M) - F Z=3 (M) - F acc (M) = F acc (M) + F mcrgc+form (M). (19) 

Each terms is defined as below, 



(20) 



poo 

F\ oca x(M) = / M'n(M',t ) dM', 

pOO 

F Z=3 (M)= M'n(M',t z=3 )dM', (21) 
Jm 

F acc (M)= [° Mn(M,t)(M) dt, (22) 

poo pto 

F acc (M)= / n(M',t)(M') dM'dt, (23) 
Jm Jt z=3 

poo pt() 

i 71 mergc+form(M) = / / M' (7mcrgc(M', t) + 7form(M', t)) dM'dt. (24) 

Jm Jt z=3 

We note that only SMBH-SMBH mergers, in which the mass of two SMBHs, Mi and M2, 
satisfy the condition that M\> M and M 2 > M, do not contribute to -F m crgc+form(M). This 
is because these mergers do not change the mass density contained in SMBHs more massive 
than M. As discussed in the previous section, -F m erge+form(-M) is not negative, thus we obtain 
the inequality, 

F Z=3 (M) < F local (M) - F acc (M) - F acc (M). (25) 

The physical meaning of each term is evident. The cumulative mass density of SMBHs 
at z — 0, F local (M), is calculated with the local mass function given by AR02. The accreted 
mass into SMBHs more massive than M during < z < 3, F acc (M), is calculated by 

F„ W - r ^ f d -^A ( ^^ *, (26) 

Jl(M) Jt(z=3) aL ec 

which resembles G acc (M, t ) defined by eq.(12). Furthermore, the mass "advection" of 
SMBHs across the boundary mass M (already mentioned in §2.1) by the mass accretion 
during < z < 3, F acc (M), is calculated by 

g,w = (1 -' )J r w r < 27 > 

ec 2 Jt(z=3) dL dM 

With the inequality (25), we can estimate the upper limit of the cumulative mass density of 
SMBHs at z = 3 (F 2=3 (M)), F local (M) - F acc (M) - F acc (M). 
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3.2. Tentative SMBHs in LBGs 

Here, we conjecture that tentative SMBHs exist in LBGs. In the local universe, the BH 
mass is not correlated with the disk mass, but with the bulge mass of the galaxy. We note 
that we tentatively treat the SMBHs in LBGs, which has no clear bulge and disk components 
(or no clear distinction between these components). Therefore, we simply assume that the 
mass of SMBHs is related to the stellar mass of LBGs as M = £M*, (see Steidel et al. 
2002). We treat the mass ratio, £, as a free parameter. To compare the proportional 
relation for the normal galaxies at z — 0, M/Mb u i ge ~ 0.001 — 0.006, we examine three cases, 
f = 0.002, 0.005, and 0.01. For the LBGs which have a typical stellar mass (~ 1O 1O M ), the 
mass of corresponding tentative SMBHs is 10 7 — 1O 8 M . The cumulative mass density of 
this tentative SMBHs in LBGs at z — 3, Flbg(^), is given by 

WM) = I , lbo( m')^-^ ? m; (28) 

where d<p(My) / dMy is the LFs of LBGs at z = 3 in the rest-frame optical (V-band) luminos- 
ity and ?7lbg(^0 is the probability that one LBG possesses a SMBH of mass M. The LFs of 
LBGs are typical Schechter-type functions and we use the one given by Shapley et al. (2001). 
Shapley et al. also analyze about 100 LBGs with a population synthesis model, assuming a 
constant star formation rate. They show the relation between the stellar mass of LBGs and 
apparent luminosity in V-band (see their Fig. 13), which is approximately expressed as 

6 A ( M* , 2 



M v = -- ^hg^h 2 ) +8.0 ) +5.01og(/i). (29) 



Here, we define the mass average of 77lbg(-^0 as 

Flbg(M) F L bg(M) 

dM 



n (M)= j lbgW _ r LBG ym) 

VhBGi } ~ r dm^dMyl - f lbg (M) m 



, M dM v dM* £ 

We reliably set F Z= ^(M) > Flbg(^), since there may be other populations of galaxies at 
z ~ 3 besides LBGs and these population may also possess the SMBHs. Very recently, Franx 
et al.(2003) found the significant population of red galaxies in z > 2, which are chosen by 
their colors, J s — K s > 2.3. They point that the number density of these red galaxies at 
z ~ 3 can be a half of that of LBGs. Furthermore, van Dokkum et al. (2003) argue that 
several among these high-z red galaxies display the AGN-like rest-frame UV spectra. More 
spectroscopy with more samples will make the AGN fraction of these high- z red galaxies 
clear. With eq.(25), we can calculate the upper limit of r/ LBG (M), 

- C\/r\ -^local(^f) - -Facc(^) - F acc (M) _ _ 

VlBg( M ) < £ 7—^ = ??LBG,maxW (31) 

^lbg{M) 
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We again stress that this is a rather loose upper limit due to some reasons. We simply 
subtract the accreted mass density during < z < 3 from the mass density of SMBHs found 
in the center of local galaxies and compare it with that of tentative SMBHs in LBGs. If 
we can know the merger rates or direct formation rates at each redshift, we can subtract 
F mergc+ f orm (M) and obtain a more strict upper limit. To investigate the redshift distributions 
of SMBH-SMBH merger rates and direct formation rates of SMBHs, the gravitational wave 
(GW)) will be a powerful tool. Today, several powerful GW detectors are under prepara- 
tion, for example, Laser Interferometer Satellite Antenna (LISA) in the low frequency band 
(1CT 4 — 1CT 1 Hz) and Laser Interferometer Gravitational- Wave Observatory (LIGO) in the 
high frequency band (10 1 — 10 3 Hz). Since both the mergers between SMBHs and the grav- 
itational collapse of SMS will emit the low frequency GW (see Haehnelt 1994 for mergers, 
Saijo et al. 2002 for collapse of SMS), these events are main targets of LISA. Furthermore, 
these any cosmological events at z < 10 should be reliably detected with high signal-to-noise 
ratio with LISA. Currently LISA is planned to start observations at the earliest in 2011. 

3.3. Upper Limit of the Probability that LBGs have SMBHs 

3.3.1. Case (i): Non-negligible Contribution of Mergers and/ or Direct Formation of 

SMBHs 

We calculate %BG,max(^0; the upper limit of the mass averaged probability that one 
LBG possesses a SMBH more massive than M, by eq.(31). Of course, ^LBG,max(^0 depends 
on mass conversion efficiency, e, and Eddington ratio, /ecM- First, we assume that these 
parameters are constants. As discussed in the previous section, in this case, the contribution 
of mergers and/or direct formation of SMBHs is not negligible and become significant in the 
low-mass end. Incidentally, many theoretical models for QSO LFs at high redshifts based 
on the hierarchical structure formation adopted the constant parameters (e.g. Hosokawa 
2002, Wyithe & Loeb 2002). Of course, even if we can not neglect the mergers and/or direct 
formation, the parameters do not have to be constant. Here, we calculate ^LBG,max(-^0 
with constant parameters as one possible case with non-negligible mergers and/or direct 
formation of SMBHs. As we have shown in the previous section, we cannot choose these 
parameters freely since the inequality (9) must be satisfied. We calculate ^LBG,max(^0 w ^ n 
the parameter sets, (e, /lbg)> indicated by points 1 ~ 3 in Fig.2. Below, we choose the mass 
range of tentative SMBHs to be 1O 7 M < M < 1O 8 - 5 M to calculate ^ LB G,max( M )- For 
M < 1O 7 M , the luminosity range of LBGs corresponds to the faint end of LFs at z ~ 3, 
thus no sufficient data is available. The fitting formulae for LFs of Shapley et al. (2001) are 
valid only for M v - 51og(/i) < -20.5. Whereas, the HXLFs of AGNs given by UAOM03 
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are valid to the sufficiently faint end, \ogL x = 41.5, which corresponds to M ~ 1O 5 M if 
we assume / E dd = 1-0. For M > 10 8 - 5 M Q , conversely, it is difficult to satisfy the inequality 
eq.(9). 

We represent the calculated % B G,max( M ) for £ = 0.002, 0.005 and 0.01 in Fig.3 with 
SMBH mass function (a). In these figures dashed lines represent the mass range where the 
luminosity of the corresponding LBG is too faint, My — 5 log(/i) > —20.5. We can see that 
rJ LBGmax (M) ~ 1 is possible for £ = 0.002 with large energy conversion efficiency (e ~ 0.3) 
(points 2 and 3 in Fig. 2) ; that is, almost all LBGs can possess the SMBHs. This value, 
£ = 0.002, is the same as the ratio between the bulge mass and the SMBH mass found with 
local galaxies. However, for £ = 0.01 only a part of LBGs can possess the SMBHs ; the 
upper limits are about 30%. If we use mass function (b), calculated 77 LB g max(^) becomes 
larger only by factor 1.5 than that in Fig.3 with similar parameter sets. We can explain why 
r/ LBG max (M) increases steeply as M decreases. As mentioned in the previous section, the 
contribution of mergers and/or direct formation to grow SMBHs increases as M decreases 
with the constant parameters of mass accretion. We do not include mergers and/or direct 
formation of SMBHs (F meise+ { oim (M)) in eq.(31) since the event rates of these processes are 
highly unknown, then the constraint on ?7lbg(-^0 becomes loose as M decrease. Therefore, 
if we can know the redshift distribution of merger rates and/or direct formation rates of 
SMBHs, we can make the upper limit more strict. Especially, if the mergers and/or direct 
formation occur mainly at z < 3, the constraint will be as strict as that with point 1 in 
Fig. 2. In Fig.3, we also plot Fi ocal (M) / F LB g(M) for comparison. This is the upper limit for 
which the population of tentative SMBHs do not exceed the local population of SMBHs. 

Fig.4 represents the mass dependence of F local (M), F LBG (M) and F local (M) — F acc (M) — 
F acc (M) respectively. These are all cumulative mass densities; F local (M) is the observational 
estimate at z — 0, F LBG (M) is one assuming that all LBGs have SMBHs with mass ratio, 
£, and F local (M) — F acc (M) — F acc {M) is the upper limit at z — 3. In this figure, we use 
the mass function (a) and adopt the parameter set displayed with the points 1-3 in Fig.2. 
It is clear that if all LBGs have SMBHs with large mass ratio as £ = 0.01, the resultant 
mass function exceed even the local estimate (F\ oca \(M)). The ratio between F^q(M) and 
F\ oca \(M) - F acc (M) - F acc (M) is our calculating ^ LBG max (M). 

We must note that we can only set the upper limit of the probability, in principle, that 
LBGs have SMBHs. Then, we can not even rule out that any normal LBGs (not AGNs) do 
not have SMBHs. Here, we conclude as follows. If we assume the mass ratio, £(= Mbh/M*), 
of tentative SMBHs in LBGs is comparable to be observed MeH/^buigc of normal galaxies at 
z = 0, we cannot completely rule out that the large part of LBGs possess inactive SMBHs. 
However, some conditions are needed in this case. That is, large energy conversion efficiency 
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(e ~ 0.3) and significant mergers and/or direct formation of SMBHs at high redshifts (z > 3). 
Whereas, if £ is larger than the local value of M-g^/ M hn \ gei the fraction of LBGs which can 
possess SMBHs decrease to less than several tens of %. 



3.3.2. Case (ii): Dominant Contribution by Mass Accretion to Grow SMBHs 

Some recent works suggest that the mass accretion may be the dominant process to 
grow SMBHs ( Yu & Lu 2004, Marconi et al. 2004). Here, we consider the special case, 
where the mergers and direct formation of SMBHs through gravitational collapse of the SMS 
are not important for M > 10 7 M Q . Concretely, G merge+ f orm (M, t ) can be neglected and the 
equality is satisfied in eq.(9), 

Giocal(M,t )=Gacc(M,to). (32) 

In this case, we differentiate eq.(32) by M and we can represent the energy conversion 
efficiency as the function of M, 

L ho \(M) [ t0 d<$> dL 



;(M) = —± J^Jt (33) 



pea 

/ n(M',t )dM' 



(YT02). Inversely, if one calculate G\ oca i(M, t ) and G acc (M, t ) using e(M) with eq.(33), 
eq.(32) is naturally satisfied in all mass range. The functional form of e(M) is the same as 
that in Fig. 4. of YT02. Here, we calculate ^LBG,max(^0 with the mass dependent energy 
conversion efficiency, e(M) given by eq.(33). We note that when mergers and direct forma- 
tion of SMBHs are negligible, the equality is satisfied in eq.(25). However, the inequality 
F Z=3 (M) > F LBG (M) still remains due to possible other populations of hosts of SMBHs 
besides LBGs. Therefore, ^lbg max(^) calculated with e(M) still gives the upper limit of 
fj LBG (M) though it is more stringent than the case with constant parameters. 

Fig. 5 represents the results. As expected, the calculated ?7lbg max(-^0 i s m ore strict 
than the results with constant parameters (Fig. 3). Even if the value of £ is comparable to 
MeH/^buigc (0.001-0.006) of local galaxies, only a part of LBGs are allowed to have SMBHs 
(less than 10-30%). Especially, if £ is larger than MBH/^buigcj LBGs which possess SMBHs 
become very rare (less than 10 %). If this is the case, 3% of LBGs which show the AGN 
activity may be the special LBGs where SMBHs are born within the galaxies. In Fig. 5, we 
also plot F local (M)/F L BG(^) for reference. These lines are same as Fig. 3. In this case, the 
difference between calculated %BG,max(^) an d F\ ocal (M) / F LBG (M) is large. In Fig. 4, we 
can see that F local (M) — F acc (M) — F acc (M) is smaller than F local (M) by about one order of 
magnitude in 1O 7 M < M < 1O 8 - 5 M . 
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3.3.3. Necessary Condition to Neglect Mergers and/ or Direct Formation 

As we show above, how much mergers and direct formation of SMBHs occur until z = 
is important to estimate the mass function of SMBHs at high redshifts (z ~ 3 in this paper). 
As mentioned above, we can observe these events by the gravitational radiation with LISA 
in the future. Here, we crudely estimate the necessary condition to drop the merger term in 
eq.(5). In eq.(2) f2(Mi,M 2 ) is the number of merger events among the SMBHs of mass Mi 
and M 2 per unit volume per unit cosmological time. We can express this f2(Mi, M 2 ) as, 

= SSi2^W (34) 

With this expression, the contribution of mergers to G mcTgc+ { OTm (M,to) may be written as 

G m ,m t „) S /J f° 1(M, M>) lmm (M>, t) dM'dt ~ ^-J^M__ . Az (35) 

Here, we transform the time interval, At to the redshift interval Az. If G merge (M, to) is less 
than G local (M, t ) and G acc (M, t ), say, at M ~ 1O 7 M , we have 

A1 ^f^,. TT7 • Az < l(T 2 (Mpc- 3 ), (36) 
A In MA In M'AzAV v ' 

where we use Gi oca i ~ G acc ~ 1O 5 M • Mpc~ 3 at M ~ 1O 7 M (see Fig.l) and we set M' ~ 
1O 7 M . Of course, other mergers than M ~ M' ~ 1O 7 M contribute to G mcrge (lO 7 M , t ), 
but we estimate loose upper limit of the event rate of mergers where M ~ M' ~ 1O 7 M to 
drop the merger term in eq.(5). We can convert eq.(36) to the number of the events which 
we can observe per unit time from the earth, 

dlnMdhM'dzdt ' dZ = A In MA\n M'AzAV ' ^ X ^ + z) ' ^ (3?) 
where d^(z) is the angular diameter distance. The time interval dt on the left-hand side is 
related to the volume interval Ay on the right-hand side as AV = 4nc(l + z)d\dt. If we 
set z = 3, the angular diameter distance is d^(z = 3) ~ lGpc in the ACDM universe (the 
redshift dependence of is not large), then we have 

dN 

— n 7T7-; , • dz < 0.6 (times • yr^ 1 ) (38) 
dlnMdlnM'dzdt j ) \ ) 

with eq.(36). This is a necessary condition for event rate of mergers of M ~ M' ~ 1O 7 M 
to drop G me rge(lO 7 M , to) in eq.(5). That is, if we observe one or more events of the SMBH- 
SMBH merger between 1O 7 M SMBHs per one year, the effect of mergers to grow SMBHs is 
at least not negligible, compared with the mass accretion. To neglect the mergers in terms of 
the continuity equation of SMBHs, no detection of events for several years is at least needed. 
We can determine the contribution of mergers observationally in the future. 
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4. Discussions 

4.1. Scatter in the M — a Relation 

In §2.2, we have searched the parameter sets, (e, /ecm) which satisfy inequality (9) with 
the HXLFs of AGNs (Fig. 2). In fig. 2, we can find the allowed region of the parameter space 
where ineq.(9) is satisfied only in M < 10 8 ' 5 M Q . However, if we demand sub-Eddington 
case, /Edd < 1-0, it is still difficult to satisfy ineq.(9) in the high-mass end (M > 1O 9 M ). 
For example, even if we adopt the most adequate parameter set, (e, /Edd) — (0.3,1.0) to 
suppress G acc (M, to), ineq.(9) is satisfied only in M < 1O 9 M (see Fig.l). One possible 
resolution of this discrepancy is the scatter included in the mass function of the SMBHs. 
That is, the scatter in the M — a relation. Fortunately, Yu & Lu (2004) study the effect of 
the scatter on the mass function. They assume the Gaussian scatter for the distribution of 
logM B H at a given a (see their eq.(43)), and find that only the high-mass end of the mass 
function (> 3 x 1O 8 M ) is sensitive to the scatter. Furthermore, the scatter (they test the 3 
cases for the intrinsic standard deviation, Ai ogMBH = 0, 0.27 and 0.4) increase the high-mass 
end (see their Fig.l). Especially, they show that the effect of the scatter is significant in 
M > 1O 9 M . This is an important result. The discrepancy discussed above will be resolved 
with the accurate decision of the intrinsic scatter in the M — a relation. Incidentally, we note 
that the fitting function of the mass function by AR02 used in this paper becomes worse in 
the high-mass end (M > 1O 8 - 5 M ) (see their Fig.8 and 9). 

4.2. Relaxing Assumptions 

In this paper, we use some assumptions for simplicity. Here, we discuss the effects of 
relaxing these assumptions. Some of these assumptions are listed as following, 

(i) extrapolation of HXLFs of AGNs to z > 3, 

(ii) no redshift evolution of e and /Edd, 

(iii) mass conservation during mergers. 

Throughout this paper, we extrapolate the observed HXLFs which are valid only at 
z < 3 to the higher redshifts (assumption (i)). However, we note that the calculations in 
this paper depend only on the HXLFs just before z — 3. To confirm this, we perform the 
same calculations with the tentative HXLFs, d&(Lx, z)/dlogL x = for z > 3.5. First, the 
critical value of e for a given / E dd to satisfy ineq.(9) becomes less than 15% smaller than that 
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of Fig. 2. Next, ?7lbg max(-^0 becomes less than 5% smaller than that of Fig. 3. Therefore, we 
conclude that the results in this paper are robust unless the HXLFs change suddenly from 
the extrapolation of the observed ones at z < 3 just before z — 3. 

Since the redshift evolution of the parameters, e and /sdd is uncertain, we adopt the 
assumption (ii). Due to this assumption, we avoid to introduce extra parameters to represent 
the redshift evolution. For the energy conversion efficiency, e, we have no idea on the redshift 
evolution. Many theoretical works adopt e with no redshift evolution (e.g. Kauffmann & 
Haehnelt 2000, Wyithe & Loeb 2002, Volonteri, Haardt & Madau 2003). For Eddington 
ratio, /Edd) several theoretical works argue that f-^dd decreases as redshift decreases and that 
this is one possible origin of the rapid decline in population of bright QSOs from z ~ 2.5 to 
z = (e.g. Kauffmann & Haehnelt 2000, Haiman & Menou 2000). Therefore, we test the 
following simple decline of /Edd to relax assumption (ii), 



We calculate G acc (M, t ) for /Edd,o = 0.5, 0.1, 0.05 and 0.01. However, the calculated G acc (M, t ) 
with /Edd(-z) is different from that with constant /sdd only by less than 10%. In this case, 
the best-fit values of constant /Edd are /sdd = 0.8, 0.64, 0.62 (corresponding to fEdd(z = 1.5)) 
and 0.5 (/Edd(^ = 1-25)) for /Edd,o = 0.5,0.1,0.05 and 0.01 respectively. 

During the merger of the SMBHs, the strong GW is emitted. Actually, the total mass of 
the SMBH binary can generally decrease due to the GW emission during the merger despite 
the assumption (iii). For the most efficient radiative merging, the total entropy (or the area) 
of SMBHs is conserved. In this case, the resultant mass of the merger between SMBHs of Mi 
and M 2 is M = y ' M\ + Mf . That is, about 30% of M x + M 2 can be maximally radiated as 
the GW. In this paper, we constrain the upper bound of the population of SMBHs at z = 3 
based on the assumption that the mass density of SMBHs never decreases. If significant 
mergers decrease the mass of SMBHs, we may not rule out that there are already much 
SMBHs at z = 3 than that at z = 0. YT02 consider the growth history of SMBHs for the 
maximally efficient radiative merging using the continuity equation. Here, we only note that 
the formulation in this paper can be applied for this case by replacing the variable, M with 
the entropy (or area) of the SMBHs, though there is the uncertainty on the spin parameter, 
a. 




for z > 2.5 
for z < 2.5 



(39) 
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4.3. Other Possible Constraint 

The current concrete observational information on the "active" LBGs is only the number 
fraction to the all LBGs, 3% ( Steidel et al. 2002, Nandra et al. 2002). However, we can 
possibly use this fraction, 3% as the lower limit of ?7lbg max(-^0- Though 3% is much smaller 
than the calculated values of ?7lbg max(^0 i n m &ny cases, if this fraction increases in the 
future, we can get the more strict lower limit. For the dominant mass accretion case (§3.3.2 
or Fig. 5), especially, calculated 77 LB g max(^) * s comparatively small, then the lower limit 
will be valid. If we can put the meaningful lower limit as noted above, we can constrain 
the lower bound of Eddington ratio, f^dd and the upper bound of the mass ratio, £. In all 
cases, calculated J? LB GminM decreases to the high- mass end. However, we must consider 
the scatter of the M — a relation in the high- mass end to apply the lower limit (§4.1). 

5. Conclusions 

In this paper, we study demographics of SMBHs in the local and high- z universe and we 
crudely constrain the possibility that LBGs have SMBHs for an application. To constrain 
the possible population of SMBHs within the LBGs, we estimate the upper limit of mass 
density of SMBHs at z ~ 3. In other words, if there are too many and too massive SMBHs 
already at z ~ 3, the resultant mass density of SMBHs at z ~ exceed the local mass 
density of SMBHs. We use the continuity equation of SMBHs, eq.(l) and transform it to the 
useful integrated form. We consider the three physical processes to grow the SMBHs; mass 
accretion, mergers and direct formation of SMBHs. Of these processes to grow SMBHs, we 
can use the HXLFs of AGNs to estimate the cosmological accretion history into the SMBHs. 

First, we investigate the demographics of SMBHs in the local universe. YT02 argue that 
inequality (9) is not satisfied with the accreted mass at z > estimated with OLFs of QSOs 
and with the local mass function of SMBHs with typical parameter set, (e, /Edd) = (0.1, 1.0). 
We use the HXLFs of AGNs including optically absorbed AGNs, then the accreted mass 
with HXLFs should have become larger than that with OLFs. However, we show that 
the discrepancy is resolved using HXLFs with the same parameters and constant Lx/^boi- 
The origin of the discrepancy is probably the assumption of the constant L-q/L^ q \. The 
comparison between the redshift evolution of OLFs and that of HXLFs suggests that Lb oc 
L]^. Therefore, if the assumption of the constant Lx/L ho \ is reasonable, the ratio of B-band 
luminosity to the bolometric luminosity has the dependence as Lb/L^i oc L^ 4 oc M~ oa . 
This makes the redshift evolution of OLFs as slow as that of HXLFs as z increases, and 
calculated G acc (M, t ) with OLFs becomes small. Furthermore, we constrain the energy 
conversion efficiency, e and Eddington rate, f-^dd to satisfy the inequality (9) upto the high 
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mass end (M < 10 8 M o , 1O 8 5 M ). If we demand /Edd < 1, the allowed region of (e, /Edd) 
exist only if we use HXLFs (e ~ 0.1 — 0.3 and / E dd ~ 0.3 — 1.0). Totally, we conclude on the 
demographics of SMBHs in the local universe as follows. 

• The cosmological accretion history based on HXLFs is more plausible than that based 
on OLFs. This is because we can satisfy ineq.(9) with typical parameters (especially 
/Edd < 1)- The reason why calculated G acc (M, t ) with HXLFs is smaller than that 
with OLFs is that we assume constant luminosity ratio, Cb and Cx for both bands. 

With these constant parameter sets, Gi oca i(M, t ) — C acc (M, t ) becomes large at the low 
mass end. This means that the contribution of mergers and/or direct formation of SMBHs 
is important in the low mass end, if the physical parameters are actually constant as we 
assumed. We adopt the constant parameters to estimate the accreted mass to SMBHs 
during < z < 3 for the case that there is non-negligible mergers and/or direct formation 
in the latter-half of this paper. 

Next, we derive the inequality (25) between the population of SMBHs at z = and z = 3 
and the accreted mass into SMBHs during < z < 3 from eq.(l) again. This gives the upper 
limit of the cumulative mass density of SMBHs at any high redshifts. For an application, we 
consider the tentative SMBHs within the LBGs. We assume that the mass of SMBHs, M is 
proportional to the stellar mass of LBGs, with the parameter £ = M/M*. We define the 
mass averaged probability that one LBG harbors a SMBH more massive than M, fj LBG (M), 
and calculate the upper limit of fJ LBG (M), ^LBG,max ( -^0 • We simply estimate the upper 
limit of cumulative mass density of SMBHs at z ~ 3 considering only mass accretion during 
< z < 3. For the event rate of mergers and/or direct formation of SMBHs, we consider 
two simple limit; (i) there are non-negligible contribution of mergers and/or direct formation 
(ii) the mass accretion is the dominant process to grow the SMBHs. The constraint on the 
population of the tentative SMBHs in LBGs are summarized as follows. 

• In case (i), we adopt the constant parameters for one possible case as mentioned above. 
If the assumed mass ratio, £ is comparable to M/M bulge observed in the local galaxies 
(£ = 0.002), we cannot rule out that almost all LBGs can have SMBHs (?/LBG,max(^0 ~ 
1.0). However, if this is the case, the large energy conversion efficiency, e ~ 0.3 and 
significant mergers and/or direct formation of SMBHs at z > 3 will be needed. If 
£ is larger than the local value (£ > 0.005), only a part of LBGs can have SMBHs 
0?LBG,maxM ~ several ten %). 

• In case (ii), we can make the upper limit more strict. In this case, we can express the 
energy conversion efficiency as the function of the SMBH mass, e(M). With this e(M), 
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calculated ?7LBG,max(^0 is onr y ten % even ^ £ * s small. Especially, if £ is larger than 
the M/Mbuige of local galaxies, ^LBG,max(^0 can become less than 10%. 

We roughly estimate the necessary condition to neglect the merger term in eq.(5). For 
example, for the major merger between SMBHs of 10 7 M Q , the loose upper limit is about 1 
event per one year. If we observe these merger events one or more times per year, we can 
not neglect the SMBH-SMBH mergers. Totally, the possibility that nearly all LBGs have 
SMBHs with large mass ratio, £ > 0.005, is reliably ruled out. Even if AGN-like LBGs 
have SMBHs with large mass ratio as Vestergaard (2003) and Vestergaard (2004) suggested, 
remaining normal LBGs will not have inactive SMBHs with the same mass ratio. 

Though we consider the tentative SMBHs in LBGs with mass ratio, £ = M/M*, this 
proportional relation is only a trial one. Papovich et al.(2001) note that the stellar mass 
estimated with the rest-frame UV-to-optical spectrum is only the mass of young stars. That 
is, additional population of old stars may exist. Papovich et al. estimate that the upper limit 
of the mass of old stars is 3-8 times as large as that of young stars. Since the SMBH mass 
is not correlated with the younger disk component but older bulge component in the nearby 
galaxies, it may be better to relate the tentative SMBH mass with the mass of old stars if 
they exist. However, our approach will be still valid, since we simply try another values of 
£. Of course, even if SMBHs actually exist in the high-z galaxies, the proportional relation 
between M and M bulge itself is never guaranteed. Fortunately, it is known that M/M bulgc of 
nearby AGNs is the same as that of local galaxies (Wandel 2002) . It will be the first step to 
investigate the relation between the stellar mass and SMBH mass of LBGs which show the 
AGN activity. 

I am grateful to S.Mineshige, Y. Ueda, T.Miyaji, A.Inoue and N.Kawakatu for helpful 
comments and their encouragement. I thank M. Vestergaard for motivating me. This work is 
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log(M BH /M ) 

Fig. 1. — Calculated G acc (M, t ) and G\ oca x(M, t ) using the local mass function of SMBHs 
(a) (see text in detail). Upper panel: G acc (M, t ) with the typical parameter set, (e, /Edd) — 
(0.1,1.0) using OLFs of QSOs and HXLFs of AGNs respectively, and Gi oca i{M,t ). Lower 
panel: G acc (M, t ) with the parameter sets indicated by points 2 ((e, /Edd) = (0.3, 1.0)), 3 
((0.3,0.3)) in Fig.2. Gi ocal (M, t ) is the same as upper panel. 
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Fig. 2. — The region of the parameter space where the inequality, Cn oca i(M, t ) > G acc (M, t ) 
is satisfied only at M < M . We calculate G\ oca \(M,to) using the local mass function of 
SMBHs (a) (see text in detail). The allowed region for each Mq is above the labeled lines 
for Mo = 10 8 ' and 1O 8 ' 5 M . We set the upper limit of the energy conversion efficiency as 
e = 0.3. The representative parameter sets are indicated by points 1 ~ 3. We can show the 
corresponding G acc (M, t ) for 1 ~ 3 in Fig.l. 
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log(M BH /M ) log(M BH /M ) log(M BH /M ) 

Fig. 3. — Mass averaged probability that one LBG possess a SMBHs more massive than M, 
r/ LBG max (M) for £ = 0.002, 0.005 and 0.01. The lines represent ^LBG,max(^0 calculated with 
parameter sets indicated in Fig.2 by points 1 ~ 3 (left to right; point 1, 3, 2 in turn). In the 
figure, dashed lines represents the mass range where the luminosity of the corresponding LBG 
is too faint : M v - 51og(ti) > -20.5. The thin solid lines represents Fi ocal (M)/F LB G(^0- 
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Fig. 4. — The mass dependence of the cumulative mass densities of SMBHs; F local (M): 
observational estimate at z — 0, F LBG (M): assuming that all LBGs have SMBHs with mass 
ratio, f = 0.01 and 0.002, at z = 3, F local (M) - F acc (M) - F acc (M): the upper limit at z = 3. 
We represent _Fi oca i(M) — F acc (M) — F acc (M) in five cases; (a): with constant parameters 
indicated by point 2 in Fig. 2 , (b): same as (a) but for point 3, (c): same as (a) but for point 
1, (d): with mass dependent energy conversion efficiency (eq.(33)) and / Edd = 1.0, (e): same 
as (d) but for f^dd = 0.3. 
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Fig. 5. — Mass averaged probability that one LBG possess a SMBHs more massive than 
M, 

iLBGmax(^) with the mass dependent energy conversion efficiency given by eq.(33) for 
£ = 0.002, 0.005 and 0.01 respectively. Dashed lines and dotted lines represents the results 
with / Edd = 1.0 and 0.3 respectively. We calculate % B G,max( M ) onl y at M = 10 7 , 10 7 - 5 , 10 s 
and 1O 8 - 5 M and connect these values with lines. The thin black solid lines represents 
F local (M)/F LBG (M). 



